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SUMMARY 


The methods of NACA Reports 865 and IO 9 O have been applied to the 
calculation of the rolling- and yawing-moment coefficients due to rolling 
for unswept wings with or without flaps or ailerons. The methods are 
based on lifting-line theory end allow the use of nonlinear section lift 
data. The method presented herein permits calculations to be made some- 
what beyond maximum lift for wings having no twist or continuous twist 
and employing airfoil sections which do not display laxge discontinuities 
in the lift curves. Calculations can be made up to maximimi lift for wings 
with discontinuous twist such as that produced by partial-span flaps of 
ailerons, or both. Two calculated examples are presented in sinqplified 
computing forms in order to illustrate the procedures involved. 


INTRODUCTION 


The calculation of the rolling- and yawing -moment coefficients due 
to rolling has received extensive treatment in the linear lift range; 
several of the better -known sources are references 1 to 3* Methods for 
making the calculations in the nonlinear lift range, however, are com- 
paratively few (for example, refs. U to 6) and are based on the over -all 
lift and drag of a wing. 

Methods for calculating wing characteristics in the nonlinear lift 
range by using lifting- line theory and nonlinear section data (refs. 7 
and 8) have given estimates which agree closely with results of wind- 
tunnel tests of nonrolling wings, as seen in references 8 and 9* It was 
believed, therefore, that these methods could be utilized in the calcu- 
lation of the rolling derivatives and that such calculations mi^t be 
more accurate than those made by existing methods. 

Although the application of the lifting-line method to a rolling 
wing in the nonlinear range is implicitly contained in reference 7 and 
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partially illustrated in reference 8, it is the purpose of the present 
report to outline the procedure of calculation by means of several 
illustrative examples. In addition, some new considerations regarding 
the application of these methods to the calculation of the rolling char- 
acteristics eire presented and discussed. Because of the supplementary 
nat\are of this report, the reader should be reasonably familiar with 
references 7 and 8. 


SYMBOLS 


■'D 


The term "section" as used herein denotes the section character- 
istics in three-dimensional flow. 


aspect ratio 

wing drag coefficient 

wing rolling-moment coefficient 

coefficient of wing damping in roll, 

wing lift coefficient 

wing yawing-moment coefficient 

coefficient of wing yawing moment due to rolling, 
effective edge-velocity factor for symmetrical part of lift 


Cn 

C 


Hp 


E 


E' 


F 

K, Ki 


distribution 




effective edge-velocity factor for antisymmetrical part of 
lift distribution. 




factor used in altering two-dimensional lift curves 

coefficients used in obtaining succeeding approximations of 
lift distribution 


Reynolds nxmiber 
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V 




c 

c 


s 


‘'do 


velocity 

two-dimensional linear lift-cvirve slope 
wing span 

local chord of wing 
root chord 

meein geometric chord, b/A 
section profile-drag coefficient 


(Cd^) two-dimensional profile-drag coefficient 

Cj section lift coefficient 

c^ai additional section lift coefficient for = 1.0 


Cj^ two-dimensional lift coefficient 


Cjmajf maximum section lift coefficient 


/c^ \ maximum two-dimensional lift coefficient 

V ^max/^ 

Cj(g) section lift coefficient for a particular flap deflection 


P 

r 


t 

y 

y* 

a 


angular rolling velocity, radians/sec 

number of intervals between points of calculation along span 
maximum thickness of wing section 
spanwise coordinate 

spanwise coordinate at end of flap or aileron 
angle of attack, deg 

correction for induced angle of attack, deg 
effective angle of attack, deg 
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^mk 

6 


e 

? 


X 

% 


induced angle of attack, deg 
angle of attack for zero lift, deg 


angle of attack for two-dimensional lift curves, deg 
angle of attack of root section, deg 
uncorrected induced angle of attack, deg 


induced angle -of -attack multiplier for asymmetrical 
distribution 


magnitude of discontinuity in absolute and induced angles of 
attack, deg 

geometric angle of twist, negative if washout, deg 
approximate average angle of twist. 



t . 1.273/1 - (1/ 

2A + 3.6 



deg 


angle of twist due to rolling motion, deg 


angle of twist at wing tip, deg 

taper ratio (ratio of tip chord to root chord) 

ratio of actixal two-dimensional lift-curve slope to theo- 
retical value of n ^/ 90 

area multiplier for asymmetrical distributions 


Vj^ interpolation multiplier 

ffjjj moment multiplier for asymmetrical distributions 

Superscript: 

* 


denotes value at end of flap or aileron 
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CALCULATION PROCEDURE 


Inasmuch as complete theoretical developments of the method used 
herein are given in references 7 and 8, only the items pertinent to the 
purpose of the present report will be given. The method is based on 
lifting-line theory and uses the effective edge -velocity factors E 
and E' to correct for the effects of aspect ratio, following the con- 
cepts of reference 10. Since E and E' as used herein are in them- 
selves only approximations for elliptic wings, it should be realized 
that the results given by the present method at zero lift will not agree 
exactly with the values which would be obtained by iising lifting-surface 
solutions for a given two-dimensional lift-ctirve slope. However, such 
differences are expected to be insignificant. 

The calculations for wings having no twist or continuous twist are 
discussed in the following section and are illustrated by means of 
example calculations for an untwisted wing designated as wing A. The 
calculations for wings with discontinuous twist require a slightly dif- 
ferent procedure and therefore are discussed in an additional section. 
The latter calculations are illustrated by means of example calculations 
for a wing with partial-span flaps designated as wing B. 


Wings With No Twist or Continuous Twist 

The rolling motion of a wing produces an antisymmetric linear aero- 
dynamic twist distribution which gives rise to the rolling derivatives. 

In the method of calculation used herein, the aerodynamic twist produced 
by the rolling motion is treated as a wing twist for pvirposes of obtaining 
the distribution of Components of the inclined force vectors are 

then found, from which the rolling derivatives are determined i 

Use of section data .- The two-dimensional data at the Reynolds n\mi- 
ber appropriate to each wing station are simply plotted against an 
effective angle, = a^E. Either the section lift coefficients c^ 

or the load coefficients c^c^b can be used, althou^ the use of the 

latter is believed to result in considerable savings in computing time. 

The load coefficients were used in the example calculations for wing A 
and are presented in figure 1. A single drag curve was used inasmuch 
as no accoxmt was taken of the spanwise variation in Reynolds number due 
to taper; the curve is presented in figure 2. The cxirves in figures 1 
and 2 are based on data taken from reference 11 for the NACA 65-006 sec- 
tion at a Reynolds number of 3*0 x 10^; however, the drag data were 
extrapolated to higher angles of attack on the basis of resiilts in refer- 
ence 12 for the NACA 64A006 section. 
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Deter»n'i nation of the lift distribution . - The effective angle of 
attack at any station y is given by 


where 


and 



The evaluation of equation (l) by the method of reference 7 is one of 
successive approximations, in which a distribution of c^c^b is 

initially assumed. From this distribution and with the values of the 
multipliers given either in table I for r = 10 or in refer- 
ence 7 for r = 20, a distribution of is obtained by using equa- 

tion (4); equation (5) can then be evaluated. With all the components 
on the ri^t-hand side of equation (2) determined, it is possible to 

find Og. Values of Cj corresponding to the values of ag are found 

from appropriate section data, after which a check distribution of 
cjc/b can be obtained. If these check values of C|c/b do not agree 
with those initially assumed, new values are assumed and the process is 
repeated until agreement is obtained. In the following examples, methods 
are indicated by which the differences between the check and initially 
assumed values are utilized to obtain the succeeding assumed values so 
that agreement is reached in a minimum ntmiber of approximations. Equa- 
tion (5) corrects the effective angle to account for the E' factor 

that is used with the antisymmetric component of loading when the sec- 

tion data are plotted against <Xg = cIqE. 
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Determination of the rolling- and yaving -moment coefficients . - 
After the distributions of Cjc/b and c^^c/b are determined ^distri- 
butions of c^^c^b are obtained by using the plots of section dataj, 

components of these distributions along the x- and z-axes in the wind- 
axes system give the following expressions: 


CxC 

b 




cLj^ c 7 c 

b V P i; b 


C7C 

— ('p 



sin(6p 

sin(6p 




( 6 ) 

(7) 


For the angles usually encountered, 

cos(ep - ^ 1 

and 


sln(,j, - oi) * . aj) 

SO that equations (6) and (7) become 


c, c 

^ 211 JLfe 

b b " b 180 VP 




(8)^ 

(9) 


It is easily seen that, for 6p = 0°, equations (8) and (9) reduce to 

the usiial expressions for the section contributions to the lift and drag 

of a nonrolling wing inasmuch as the contribution of Ch to the lift in 

°-o 

equation (9) is usually negligible. 

Spanwise integrations of equations (8) and (9) yield the yawing- 
and rolling-moment coefficients as follows: 
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£ 

b b b 


( 10 ) 





[' 


I + 



£ d^ 
b b b 


( 11 ) 






Jt 

l8o 



“i)m 




(11a) 


Values of the numerical Integrating multipliers 
for r = 10 and r = 20 in addition to the area 
Tjjjj multipliers may be used to find Cj^ and Cp 

that is, Ct = A ^ f— and so forth. 
m=l \ ^ /m 


Ojjj 5ire given in table II 
multipliers The 

by similar relations, 


Application of the method . - In order to illustrate the use of the 
method, an example is presented for wing A at an angle of attack of 12°, 

rolling at such a rate that the wing-tip helix angle pb/2V is 0.01 radian. 

The calculations ane made for r = 10. The pertinent dimensional data are 
given in table III and the lift distribution is calculated in table IV. 

The initially assumed load distribution in the linear lift range can 
be most rapidly obtained by using methods with which the reader is prob- 
ably already familiar. For example, nearly exact initial assumptions can 
be made very quickly by proper use of the numerous lifting- surface solu- 
tions which are available in chart form and thus effect appreciable 
savings in computing time. For the sake of presenting a numerical pro- 
cedure, however, the method of reference 13 , which has been shown to give 

very good results in most cases, was used in this report. By combining 
the equation on page 8 of reference 7 (expressed in terms of CjC/b^ with 

equation (33) of reference 8 and equations ( 13 ) and (l8) (modified for 
E' ) of reference 13, the following expression is obtained for the ini- 
tially assumed load distribution: 
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(¥) 


assumed 


1 

— 

2 


Ac 

— + 1.273 
b 




«=l(a) Aao(e - i) 


A + 1.8 


AE + 6 


Aao^p 
AE' + 4| 


( 12 ) 

where Ci(a) section lift coefficient for the geometric angle of 

attack in question. This loading is entered in column (|) of table IV 
in normal order and in column (§) in reverse order. The loading is 
entered in this manner in order to shorten the size of the conq)uting 
form. The mechanics of conqputing are all self-explanatory and it will 
be noted that the correction to the antisymmetric part of the loading is 

made in colimm (||) as obtained by using equation (5). For the sake of 
brevity only the final calculations have been shown in the table. The 
check load coefficients in colvunn @ , read from the section plots, 
will usually not agree with the assumed values for the first approxi- 
mation. The process is repeated rnitil agreement is obtained between 
the assumed and check values of c^c^b. The manner used in determining 

the succeeding assumptions is dependent on a number of factors such as 
the linearity and the slopes of the section loading curves and on 
whether r = 10 or r = 20. Various methods for obtaining succeeding 
assumptions are presented in the appendix. When calculations are being 
made for more than one angle of attack, the first assumption for the 
second angle of attack can be based on the solution for the previous 
angle of attack by finding the value of c^c/b corresixjnding to an 

angle 




2 



'ijdag 




(13) 


doe 

in which it can be assumed that 0.7 in the linear reinge. Once 

the values of have been determined for two angles of attack, plots 

of Og against ag can be made for each section. Values of c^cyb 
corresponding to the extrapolated values of will usually give a 

fairly accurate first assumption and thus minimize the amount of com- 
puting required. If only limited calculations are being made, it is 
recommended that any calculations in the nonlinear range be based on 
the results of a calculation in the linear range in a manner similar 
to that just described. In general, this procedvire will eliminate a 
rather arduous solution since the load distributions may change very 
rapidly in the nonlinear range. 
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After the induced angles of attack and the lift distribution are 
determined, then the profile-drag distribution can be determined. 
Inasmuch as each section is assumed to be acting two -dimensionally, 
the section drag coefficients are obtained at the section lift coef- 
ficients or effective angles for the proper values of Reynolds 

and Mach n\ambers. The calculations are carried out in columns (2) 
to of table V. The calculations leading to the rolling derivatives 
are carried out in col\mins © to @ and at the bottom of the table. 


Wings With Discontinuous Twist 


The discussion in this section is limited to the case where r = 20; 
a similar method could be devised for r = 10. It is emphasized that the 
method of altering the two-dimensional data, which is subsequently 
described, applies only up to and including maximum lift and therefore 
precludes the calculation of and beyond maximum lift. 


Alteration of two-dimensional data ..- The two-dimensional data to 
be used for wings with discontinuous twist must be altered in order to 
avoid a discontinuity in the spanwise distribution of maximum lift coef- 
ficient at the end of the flap or aileron since no such discontinuity 
exists in the physical flow. The maximum lift-coefficient values are 
altered (following ref. 8) by the relation 


^imax (^imax)Q ^ max) 


(14) 


where the spanwise variation of F is given for several cases in fig- 
ure 3. An illustration for deriving the F factors for any wing is 
given in figure 4. The quantity Ac^* is the increment in c^ 

at the discontinuity due to the deflection of the flap or aileron for 
the proper local Reynolds number. The values of Cj and are then 

altered according to the eqtiations 



(15) 




= + E(a^ 



(16) 
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The data altered in this manner are shown plotted in figure 5 for 
wing B with 60-percent-span split flaps deflected 60®. For purposes 
of con^JQrison the unaltered two-dimensional section data cross -plotted 
from reference 11 are also shown in figure 5. 

The two-dimensional drag values as such are not altered, hut the 
values of corresponding to the values of c^^ are replotted 

against either c^ or ttg. The drag data used in the example for 
wing B are shown in figure 6 plotted against c^. The data for the 

unflapped sections were taken from reference 11, cross -plotted to the 
values of R shown in figure 5. Since no drag data were available for 
the NACA 64-210 section with split flaps deflected, the data for the 
NACA 23012 airfoil section (ref. 14) were used inasmuch as the lift 
curves were nearly identical to those for the NACA 64-210 section up 
to maximum lift. Since the data of reference 14 were for R = 3.5 x 10^, 
no accoxmt was taken of the Reynolds nxanber variation across the flapped 
portion of the span. The manner of alteration just described is neces- 
sarily arbitrary and further experimental work may indicate a different 
proced\ire; however, it should be recognized that the drag contributions 
depend on the differences in between the right and left wings and 

therefore are not critically dependent on the absolute values of the 
drag polars . 

Determination of the lift distribution .- For wings with discontinuous 
twist the induced angle of attack computed by means of eqiiation (4) is 
modified by a term which is proportional to the magnitude of the discon- 
tinuity and acts as a correction factor to account for the inability of 
a limited trigonometric series to represent adeqiiately the lift distri- 
bution of a wing with partial-span flaps or ailerons deflected. This 
induced angle, according to reference 8, is 

% ’ % ^ KtI 


where is the uncorrected induced angle of attack given by equa- 
tion (4), 6 is the magnitude of. the discontinuity, and values of the 

correction factors taken from reference 8 are given in 

table VI. These correction factors depend only on the spanwise posi- 
tion of the discontinuity and apply to a deflected flap or aileron 
extending from the point of discontinuity to the right wing tip. Values 
for any spaiwise extent of the deflection of the flap or aileron, or 
both, can be obtained by combinations of the values presented. For 
example, if the flap deflection is symmetrical and extends over 60 per- 
cent of the wing span, the values of would be obtained by 
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2v* 

subtracting the values of a^,/6 corresponding to = 0.6 from those 


2y* _ 


= -0.6. The tabular values of 


which correspond to 
2y* 2v 

-T— = are for the flap side of the discontinuity 

D D 


a^/5 


for 




there- 


fore the corresponding section lift curves should be used with these 
values. The value of 6 is obtained as shown in figure 5. 


If the discontinuity is located at a station other than ^ = cos — , 

b r 

the values of ^c^c^^b) are interpolated from the calculated values of 
CjC/b as follows; 



(18) 


Values of the interpolation multipliers are given in table VII for 

various spanwise positions of the end of the flap* Multipliers for 
positions of 2y*/b other than those tabulated can be calculated from 
the relations given in reference 8. 

Application of the method .- The method is applied to wing B with 
60-percent-span split flaps deflected 60° at ag = 10° and rolling at 

such a rate that the wing-tip helix angle pb/2V is 0.01 radian. The 
geometric characteristics of the wing are given in table VIII and the 
calculation of the spanwise lift distribution is given in table IX. 

The determination of the initial assumption for the lift distribution 
is not given herein inasmuch as the procedure is fully illustrated for 
this wing in reference 8 for the symmetrical case. For the asymmetrical 
case the approximate antisymmetric loading given by equation (l8) of 
reference 13 (modified for E' ) 





P 


AE' + U 


should be added to the assumed symmetrical loading given by the pro- 
cedure used in the exan5>le of reference 8. Any suitable method for 
obtaining the assumed loading may be used, however. 
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The values of in column of table IX are those for 

2 y* 

= -0.6 from table VI; whereas the values of in column 

2 v* 

are the negative of those for — ^ =0.6 ajid were obtained by the pro- 

cedure outlined in the previous section. The values of S at ±2y*/b 
are found from the flapped and unflapped section curves (fig. 5 ) at 
values of Cj* corresponding to the values of ^c^c^b^*. After these 

approximate values for 6 are determined, the check span load distri- 
bution can then be calculated. The check values of ^Cjc/b)* are then 

interpolated from equation (18). When convergence is established 
(i.e., colvunn minus column is zero), the values of at 

t2y*/b can be found. (The method given for case I of the appendix 
was used to obtain convergence between the assumed and check values of 

Cjc/b, columns (3) and , respectively. ) 


The drag distribution is calculated in columns to of 
table X. Note that two values are given at — = +O.60. The values at 

2y ^ ^ 

-g- - +0.60 + 0 correspond to the flapped sides. 


/ 


The rolling- and yawing-moment conq)onents are calculated in 

col^s © to @ and the coefficients and derivatives are calculated 
at the bottom of table X. The multipliers % used in table X are 


tabulated in table XI. If the values of and Cp are desired, the 

numerical integrating area multipliers of reference 8 may be used. 


DISCUSSION 


Lack of either experimental wing-alone rolling data or suitable 
two-dimensional section data prevents the making of exact correlations 
of this method with experiment in the vicinity of maximum lift. It was 
possible, however, to con5>are the wing-alone calculations for wing A 
with experimental wing -body results, since sufficient section data 
existed to allow the calculations to be carried slightly beyond maximum 
comparisons are presented in figure 7. Agreement is con- 
sidered to be good when the differences between the conditions of the 
calculations and tests are considered. The failure of the present 
method to predict the increase in Cj for 5° < a < 8° may be 
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inherent in the method since similar experimental trends have been 
observed for other wings of this plan form. The differences between 
the calculated and experimental variations of Cn^ in the same angle- 

of -attack range could be partly associated with the previously noted 
increase in and partly due to body-interference effects as shown 

in reference 15. The differences in the hi^ lift range, particiilarly 
with regard to the variations, probably can be all ascribed to 

body-interference effects since the differences are similar in trend to 
those of reference 15. 

Also shown in figure 7 are the variations of C, and C„ calcu- 

lated by means of references 5 and 6, respectively, in which use is made 
of the experimental wing-body lift and drag curves. It appears that both 
the present method and that of reference 5 give reasonable estimates of 
the variation of The present method, however, appears to give a 

more realistic picture of the variation of 0^^ in the high angle-of- 

attack range than that of reference 6. The edge-suction effects which 
are applied in reference 6 were negligible for this aspect ratio and 
therefore do not account for the differences shown. 

Several items of interest were observed diiring the covirse of the 
calculations for wing A. For example, the contributions of the section 
lift and profile drag to the rolling derivatives coxxld be separated as 

shown in figure 8. The profile-drag component of Cv, is seen to be 

P 

opposite and nearly two to three times that due to the lift at hi^ 
angles of attack. The profile -drag contribution to the damping in roll, 
on the other hand, is either zero or negligible even at the hi^er angles 
of attack. 

The method, presented herein is primarily intended for use in the 
nonlinear range where some flow separation is present, and the subsequent 
discussion briefly covers some considerations which limit the use of this 
method. The condition that must be essentially fiilfilled is in keeping 
with the basic assimjption that all sections operate two-dimens ionally or 
nearly so. Therefore any separation that is present must not give rise 
to large amounts of spanwise flow which can cause a complete redistri- 
bution of lift and thus invalidate the basic asstimption. lack of suf- 
ficient experimental data on load distributions beyond maximum lift 
prevents the formulation of positive limits concerning these phenomena, 
but, on the basis of airfoil characteristics and on observed stalling 
behavior, regions most likely to be amenable to calculation can at least 
be classified. 

It is believed that departures from the two-dimensional-flow assunp- 
tlon will not seriously affect the calculation of the rolling derivatives 
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for wings which incorporate airfoil sections that exhibit gradual 
c ^ges in lift-curve slope beyond maximum lift. On the other hsnfi 

tion^ section characteristics of neighboring sec- 

display a tendency toward asymmetrical 

tta mounts, 

curve peaks at naxIsuM ufr^f 

iispsss-is 

airfoli*se^orstrn‘’^°"\“^ *“ ‘‘^aslble on the basis of 

slderaU^rSlch^^f Slf r%1^^=t”o?’ tJ^:^ =- 

referred to In refeLce Se^^lbllSj iSr^s”^- 


dttj 

da 


< -1 


( 19 ) 


?kTo?th?li;" 

el«aef T!lo°"r'"r “ o”^rL"’'.f calculations with 

from ewperlenoe that calculations beyond“^lmS'‘nffwUh ”r^ So”*^ 

cSSaSonrSJfScessS) ‘^o^^ver^^r^ “■ 

from calculations with /L lo be wli as th?lnlM.f“* 

for thp r - 90 e,rai,.j-- iHitial approximations 

co;;,X time ■ S'— UJ' Shorten 


l*^^Sl®y Aeronautical Laboratory^ 

National Advisory Committee , for Aeronautics, 
Langley Field, Va., January 29, I955. 
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APPENDIX 

METHODS FOR OBTAINING SUCCEEDING APPROXIMATIONS FOR THE 
SPANWISE LOAD DISTRIBUTION FROM AN INITIALLY 
ASSUMED LOAD DISTRIBUTION 
Msthods for Use With r = 20 


These methods for obtaining succeeding approximations for the load 
distribution have been extensively used and produce convergence in 
reasonably few approximations. 

dci 

Case I: positive and linear.- For a positive and linear lift- 

da 

crirve slope, the succeeding approximation, denoted by the superscript 1, 
is given by the equation 


where 



(Al) 


(A2) 


the differences between the check 

@ minus column (^ , table IX ), 
and K and are constants for any particular wing. Values of K 

and taken from reference 8 are presented in figure 9 as functions 

of AE/q. Values of for i greater than 3 are small enou^ to be 

considered negligible. 

The number of terms of equation (A2 ) needed for any particular 
approximation depends upon the convergence of the assuniptlons; fewer 

/ClC\ 

terms are needed when the differences ^ j are small or when 


in which the increments Al 


ff) 


are 


values and the assumed values (coliomn 
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positive differences nearly cancel negative differences. Equation (A2) 
applied to wing B at m = 6 = O-SSs) becomes 




+ 0.3 



For m = 1 and 2 , equation (A2 ) is expanded as 
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For m = 18 and I9, the values of 
similar in form to the values of 



For wing B, equation (A 3 ) would he 


A 



1 

8.3 


(3.2' - o.i^)(© -(3))^ + (1.0 - 


0.3)(© -0)2 + 


0.u(@ -0^ + 0.3(@ 



whereas equation (Ak) would become 


A' 





+ 


O.k 




+ 



Use of the K factors has been found to establish convergence 
within three approximations when the initial assumption is reasonably 
close. 


dci 

Case II: positive and nonlinear.- Althou^ the factors in 

da 

figure 9 were derived for a linear lift-curve slope that is constant 
across the span, they can be used in the nonlinesir range. An estimation 
is made of the wing lift-cxirve slope at the angle of attack in question, 
from which a value of q is obtained. The K factors are then obtained 
at the proper value of AE/q. In this way it is possible to use these 
factors for values of q as low as O.3 for A = 3 * 0 * From the trends 
of the variations, it is seen that considerable extrapolation is per- 
missible at the hi^er values of AE/q. 

dcy 

Case III: negative and linear or nonlinear.- For a negative 

lift-curve slope, either linear or nonlinear, the procedure is as follows: 
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(1) First obtain values of the span load distribution with fewer 
significant figures than desired. 

(2) Find at each station, and work with these values 

directly. 

(3) Adjust the load distribution at those stations where the 
largest values of A^— occur by adjusting the loadings at m and 

m ± 1 stations to obtain approximate convergence at the mth station. 

If large values of a(-~^ occur at two adjacent stations, then adjust 

the loading at only one of the stations, although with some practice 
the adjustments required at two such stations become quite obvious. 


(4) As the values of 

significant figures in the 
amount. 


A ^-^ j are made smaller, the number of 
solution can be increased to the desired 


(5) Keep all values of either positive or negative, if 

possible, for easier mental manipulation. 


Msthods for Use With r = 10 

Although factors similar to the K's given in case I could be 
derived for the corresponding condition with r = 10, it is equally 
facile to use the method described under case III for all cases with 
r = 10. 
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TABLE I.- nroUCED-ANGLE-OF-ATTACK MULTIPLIERS FOR 

ASYMMETRICAL LIFT DISTRIBUTIONS WITH r = 10 ^ 


^ m=l' ° /m 


b 

b 

-0.9511 

-0.8090 

-0.5878 

-0.3090 

0 



m 

9 

8 

7 

6 

5 



-0.9511 

9 

231.766 

-43.879 

0 

-2.148 

0 

1 


-.8090 

8 



-34.405 

0 

-2.573 

2 

.8090 

-.5878 

7 

0 


88.527 

-29.824 

0 


.5878 

-.3090 

6 

-6.610 

0 

-35.060 

75.305 



.3090 

0 

5 

0 

- 4.377 

0 

-30.000 



0 

.3090 

k 

-1.716 

0 

-3.388 

0 



-.3090 

.5878 

3 

0 

-1.188 

0 

-2.882 

0 

7 

-.5878 

.8090 

2 

-.5ifk 

0 

-.863 

0 

-2.573 

8 

-.8090 

.9511 

1 

1 

■ 0 

-.286 

0 

-.558 

0 

9 

-.9511 



1 

2 

3 

4 

5 


L Vk 



0.9511 

0.8090 

0.5878 

0.3090 

0 

b 

D 


Values of k at top to be used with values of m at left side; 
values of k at bottom to be used with vsQ.ues of m at right side. 
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TABLE II. - WING-COEFFICIENT MULTIPLIERS FOR WINGS 
WITHOUT DISCONTINUOUS TWIST 



II 

CO 

0 

r = 10 

m 


% 

m 

Tm 

% 

-0.9877 

19 

0.01638 

-0.00809 




-.9511 

18 

.01618 

-.00769 

9 

0.06472 

-0.03078 

-.8910 

17 

.04754 

-.02118 

- 



-.8090 

16 

.03078 

-.01245 

8 

.06154 

-.02489 

-.7071 

15 

.07405 

-.02618 

- 



-.5878 

14 

.04236 

-.01245 

7 

. 16943 

-.04980 

-.U 540 

13 

.09331 

-.02118 

- 



-.3090 

12 

.04980 

-.00769 

6 

.09959 

-.01539 

-.1564 

11 

.10343 

-.00809 

- 



0 

10 

.05236 

0 

5 

.20944 

0 

.1564 

• 9 

.10343 

.00809 

- 



.3090 

8 

.04980 

.00769 

4 

.09959 

.01539 

.4540 

7 

•09331 

.02118 

- 


— - — 

.5878 

6 

.04236 

.01245 

3 

.16943 

.04980 

.7071 

5 

.07405 

.02618 

- 



.8090 

4 

.03078 

.01245 

2 

. 06154 

.02489 

.8910 

3 

.04754 

.02118 

- 



.9511 

2 

.01618 

.00769 

1 

.06472 

.03078 

.9877 

1 

.01638 

.00809 

- 






2 k 


NACA TN 2957 


TABLE III.- GEOMETRIC CHARACTERISTICS OF WING A 


Taper ratio. 
Aspect ratio. 
Root section 
Tip section 


X o.6o 

A 4.00 

NACA 65-006 
NACA 65-006 


Geometric twist, €^, deg None 

Edge-velocity factor, E l.ll8 

Edge -velocity factor, E' 1.4lU 

Wing Reynolds number, R 3*0 x 10^ 


2y 

b 

c 

b 

R 

(assumed) 

t 

c 

e 

€ 

0 

0.3127 

3.0 X 10^ 

0.06 

0 

0 

• 309 

.2740 

3.0 

.06 

0 

0 

.588 

.2392 

3.0 

.06 

0 

0 

.809 

.2115 

3-0 

.06 

0 

0 

.951 

.1937 

3.0 

.06 

0 

0 


c ^ 2 

^ A(1 + X) 

Alter values of ^ near tip if 

D 

tip is rounded.) 



1 - (1 - ^ 
D 




TABLE IV.- CALCUIATION OF LIFT DISTRIBUTION FOR WING A 


NACA TN 
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TABLE V. - CALCULATION OF THE ROLLIHG DERIVATIVES FOR WHIG A 
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TABLE VII.- INTERPOLATION MULTIPLIERS v„, FOR OBTAINING VALUES 
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TABLE IX.- CALCULATION OF UPT DISTRIBUTION FOR WING B 
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^-X@«i8rnk for 19^ kg lOj. 0 |^ - ? 0i/9mk for lOj 

«2<D*0. 

[(7.49 - 4.60) - (8.63)(6.1S)] - -0.02. 
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TABLE X.- CALCULATIOR OP OHE ROLLIHG DERIVATIVES FOR WIRO B 
|^.6-spaii spilt flaps deflacted 60°; Og = 10.00®; ^ » 0.01 radla^ 


© © 

2y 

^ (table 

” IX) 


© © 


1 o 

(« 8 . 6 ) 


0x00.573x0 


© © © 


4.50 -5.07 -0.0016 


2.21 -2.76 
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(^P “ ®l) XI) 
® X © 
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TABLE XI.- ROLLING- AND YAWIHG-MOMENT-COEFFICIENT MULTIPLIERS FOR WINGS WITH 
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Figure 1.- Section loading curves used in the example calculations 

for wing A. 
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having discontinuous twist. 












Figure 6.- Section drag curves used in the example calculations for wing B, 
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Figure 8.- The calculated rolling derivatives for wing A showing the lift 
and profile-drag components. R = 3.0 X 10^. 












